Competition between quantum-liquid and electron-solid phases in intermediate 

Landau levels 
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On the basis of energy calculations we investigate the competition between quantum-liquid and 
electron-solid phases in the Landau levels n = 1, 2, and 3 as a function of their partial filling 
factor v. Whereas the quantum-liquid phases are stable only in the vicinity of quantized values of 
v = l/(2s + 1), an electron solid in the form of a triangular lattice of clusters with a few number 
of electrons (bubble phase) is energetically favorable between these fillings. This alternation of 
electron-solid phases, which are insulating because they are pinned by the residual impurities in 
the sample, and quantum liquids displaying the fractional quantum Hall effect explains a recently 
observed reentrance of the integral quantum Hall effect in the Landau levels n — 1 and 2. Around 
half-filling of the last Landau level, a uni-directional charge density wave (stripe phase) has a lower 
energy than the bubble phase. 

PACS numbers: 73.43.Nq, 73.20.Qt 



INTRODUCTION 



Two decades after the discovery of the integral and 
fractional quantum Hall effects (IQHE and FQHE), two- 
dimensional (2D) electron systems in a perpendicular 
magnetic field remain a field of research with unforeseen 
and surprising phenomenal Experimental findings range 
from the observation of large anisotropies in the longi- 
tudinal magneto-resistance^ to an intriguing reentrant 
IQHE (RIQHE) at moderate magnetic fields^ From 
a theoretical point of view, these systems are particu- 
larly interesting because they exhibit a large variety of 
quantum phases ranging from quantum liquids respon- 
sible for the FQHE to electron-solid phases such as the 
Wigner crystal and charge density waves (CDWs) . These 
quantum phases arise because in the presence of a mag- 
netic field the electrons' kinetic energy is quantized in 
equidistant energy levels, which are called Landau levels 
(LLs). The separation between these highly degenerate 
levels is fiuc, where u>c — eB/m is the cyclotron fre- 
quency for electrons with charge e and band mass m. 
The level degeneracy is characterized by the LL filling 
factor v — n e i/riB, where n e i is the planar electronic 
density and Ub = l/(2irl B ) = eB/h is the flux density, 



with the magnetic length I b — y/h/eB. The flux density 
determines the density of states per LL. Furthermore, 
each LL is split into two spin branches with an energy 
gap Az because of the Zeeman effect. Correlation effects 
due to the Coulomb interaction between the electrons 
become important if v ^ N, with integral N, and if the 
characteristic Coulomb energy e 2 /e/s is smaller than the 
level separations Tiloq and Az, which is the case at large 
magnetic fields. 

In the two lowest LLs, strong correlations give rise 
to quantum-liquid phases, which display the FQHE at 
v = pj (2ps ± 1), where p, s are integers and v — v — N is 
the filling factor of the last, partially filled spin branch of 
the n-th LL. TV denotes the number of completely filled 



levels, which may be treated as inert, and obeys N = 2n 
for a partially filled lower spin branch and N — 2n + 1 
for the upper spin branch of the n-th LL. The FQHE 
is understood in terms of composite fermions (CFs)j^ 
which arc formed to reduce the repulsive Coulomb in- 
teraction between the electrons. This reduction is due to 
a factor ni<j( z i — z j) 2s m the iV-particle wave function, 
where Zj — Xj + ijjj is the position of the j-th particle 
in the complex plane. Furthermore, these CFs experi- 
ence a reduced coupling (eB)* = eB/(2ps ± 1) to the 
external magnetic field, and thus fill CF-LLs with a fill- 
ing factor v* = n e ih/(eB)* as do electrons in the case 
of the IQHE. For p = 1, the CF wave functions coincide 
with Laughlin's trial wave functions, which describe the 
FQHE at v = l/(2.s + l). 6 Both the IQHE and the FQHE 
are usually explained by single-particle localization due 
to underlying impuritiesi if the filling factor is slightly 
increased from v = N (y* = p), e.g. by lowering the 
magnetic field, electrons (CFs) are promoted to the next 
(CF-) LL. They first populate the minima of the under- 
lying impurity potential and are thus localized so that 
they do not contribute to the electrical transport. The 
longitudinal and the Hall resistances therefore remain at 
their original values. This gives rise to the plateaus in 
the Hall resistance at values R xy = h/e 2 N for the IQHE 
(R xy = (h/e 2 )(2ps ± l)/p for the FQHE) if plotted as a 
function of the magnetic field accompanied by zeros in 
the longitudinal magneto-resistance. 

However, the insulating behavior of electrons in the 
last LL may also be due to another effect than single- 
particle localization. Because of their mutual Coulomb 
repulsion, the electrons in the last LL have a tendency 
to form a Wigner crystal at low v. At larger values 
of P, Hartree-Fock calculations indicate that in higher 
LLs these electrons may form CDWsA2*iS There is since 
quite long experimental evidence for the formation of a 
Wigner crystal in the lowest LL below v ~ 1/5. — Recent 
microwave experiments revealed that the Wigner-crystal 
phase around v = N persists in higher LLsi£ and that a 
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triangular CDW (bubble) ground state is formed around 
v = 1/4 and 3/4 in the n = 2-LLm& Another experi- 
mental evidence for CDW formation in higher LLs is the 
huge anisotropy observed in the longitudinal magneto- 
resistance at half- filling for n > 2. - In this case, the bub- 
bles merge to form a uni-directional CDW (stripe phase) . 
Because the electronic transport is easy along the stripes 
(low-resistance direction) and strongly suppressed across 
them, this effect gives rise to the anisotropy. 

Recent experimental investigations by Eisenstein et al. 
in both spin branches of the n — 1-LL have revealed a 
new intriguing phenomenon^ between the FQHE states 
at v — 1/5,1/3 and the even-denominator stated at 
v = 1/2 with a quantized R xy = h/e 2 (N+v), the Hall re- 
sistance jumps back to values R xy — h/e 2 N, correspond- 
ing to the neighboring plateau of the usual IQHE. An 
analogous reentrant IQHE (RIQHE) had been observed 
before by Cooper et al. in n = 2 around v = 1/4 and 
3/4i£ As expected from the observation of the FQHE, 
which indicates the importance of strong correlations in 
this range of the filling factor, such insulating phases of 
electrons in the last LL are unlikely to be due to simple 
localization of single particles. 

In this article, we present detailed energy calculations 
of the different quantum-liquid and electron-solid phases 
in n = 1, 2, and 3, which allow us to identify the electron- 
solid phase in form of a triangular CDW as the origin of 
this insulating behavior. The CDW is pinned by the 
underlying impurities and therefore does not contribute 
to the electrical transport. This leads to the same dc 
response in transport measurements as observed in the 
IQHE regime. However, in the vicinity of V = l/(2s + 1) 
the quantum-liquid phases have a lower energy than the 
electron solid, provided that v is lower than a critical fill- 
ing P c (n) oc l/(2n + 1), given as a function of the LL 
index This gives rise to the observed alternation 
between the different phases with first-order quantum 
phase transitions between them. The phases in n = 1 
and n — 2 have already been discussed in a previous 
Rapid Communication. 17 Here, we provide more details 
of the calculations including stripe phases and results for 
n = 3. We further evaluate how pinning by impurities 
will affect the energy values of the electron-solid phases 
and discuss possible experimental consequences of our 
theoretical investigations. 

In Sec. II, we present a theoretical model for the de- 
scription of electrons restricted to the n-th LL and its ba- 
sic properties. The Hartree-Fock solutions of this model 
for triangular and uni-directional CDWs are obtained in 
Sec. III. The energy of the quantum-liquid phases exactly 
at v — l/(2s + 1) is calculated in Sec. IV with the help of 
sum rules imposed on Laughlin's trial wave functions. 18 
If the filling factor is moved away from precisely these 
values, quasiparticles/-holes are excited. Their energies 
are calculated analytically in Sec. V in the framework of 
the Hamiltonian theory proposed recently by Murthy and 
ShankariiS In Sec. VI, a comparison of the energies of the 
different phases and the relation to experimental results 



are given. Impurity effects are investigated in Sec. VII, 
and possible experiments, which could verify our results, 
arc discussed in Sec. VIII. A brief summary may be found 
in Sec. IX. 



II. THE MODEL 

We adopt a model of spinless electrons and restrict 
their dynamics to the n-th LL. This model is valid when 
the characteristic Coulomb-interaction energy e 2 /els is 
smaller than the LL separation fkuc and the Zeeman 
splitting Az- We further require that the partial filling 
factor of the last level v = v — [v] is different from zero 
because at integral fillings, the only possible low-energy 
excitations are intcr-LL excitations, which cost an energy 
of order hue or A z . In the absence of inter-LL excita- 
tions, the number of electrons N e i in the last LL is fixed 
by the partial filling factor N e i = DA/2irl B , where A is 
the total area of the sample. The kinetic energy is there- 
fore an unimportant constant E^in — N e iTiujc{n + 1/2), 
which may be omitted in the Hamiltonian. The model is 
thus given by 

where v(r) is the usual Coulomb potential and the 
fermion fields contain only components of the n-th LL 

^«( r ) = ^2( r \ n >yo)Cn,y - 

yo 

The wave functions in the Landau gauge A = B{—y, 0, 0) 
are given by (r\n,y ) = L" 1 / 2 exp{iy x / l 2 B )\ n {y - y ) 

with X n(y) = {^2 n n\l B )- 1 ' 2 H n {y/l B )eM-V 2 l^l): 
and c\ yo creates an electron in the state |n, yo). In terms 
of Fourier components of the density operator 

Pn(q)= / d 2 r^(r)^(r)e lq - r = F„( g )p(q), 

the Hamiltonian can be rewritten in reciprocal space as 

q 

where = A J d 2 q/(2n) 2 . The projected density oper- 
ator reads 

=Y. e ~ lqVV ° C iyo+ q *ll/2 C n,yo-< 1 *ll/2> 

yo 

and F n (q) = L n (q 2 l 2 B /2)exp(-q 2 l 2 B /4), where L n (x) is 
a Laguerre polynomial. The form factor F n (q) is due to 
the overlap of the wave functions in the n-th LL and has 
been absorbed in an effective interaction potential 

v n (q) = v(q)[F n (q)} 2 , v(q) = 2ire 2 /eq. (2) 
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FIG. 1: Effective interaction potential v(r) in real space in 
units of e 2 /els- Points correspond to n — 1, 5; gray line: 
pure Coulomb potential; black line: scaling form according to 
Eq. @. 



Inter-LL excitations may be accounted for in a wave- 
vector-dependent dielectric function e(q)&2& This depen- 
dence will be neglected in the further discussion because 
it has only a minor effect on the physical properties of 
the system as shown by Fogler and Koulakov. 10 

With the help of the usual anti-commutation relations 
for the fermionic operators, {cn, yo , 4i>, y '} = <>n,n<$y ,y' 

and {c n ,y ,c n > t y> } = {4,bo ' c n',^ } = °' one finds tnat 



the projected density operators satisfy the algebri 
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[p(q).p(k)] = 2isin 



(q x k)J% 



p(q + k). (3) 



The physical origin of this algebra is the following: by 
projecting the density operators to the n-th LL, one av- 
erages over the rapid cyclotron motion so that p(q) may 
be interpreted as a density operator of fermions, which 
are described only by their guiding-center coordinates 
R s ; = r s ; — z x Pi/eB with and being position 
and momentum of the i-th particle, respectively. It is 
the non-commutativity of the components of the guiding- 
center coordinates [Afc,Yj] = —il B 8k,i, which is respon- 
sible for these commutation relations as well as for the 
fact that each electronic state occupies a minimal surface 
a = 2nl B = l/n B . 

The Hamiltonian Q together with the commutation 
relations (PJ define the quantum mechanical model. The 
electrons interact via the effective interaction potential 
(0), which contains the information of the LL of inter- 
est in the form factor F n (q). Deeper insight into the 
meaning of this interaction potential is obtained from 
a transformation back to real space, which is shown in 
Fig-H The effective interaction potential satisfies the 
scaling law v n (r) = v(r/Rc)/\/2n + 1, where Rc = 
Ib\J2u + 1 is the cyclotron radius^ With the help of 
F n (q) ~ Jo(q^/2n + 1), which becomes exact in the large- 



n limit, one obtains 
4e 2 

v(x) ~ — ; — Re 



itcIbx 



K 



1 - yjl - A/x 2 



(4) 



where Jo (a;) is the zero-order Bessel function and K(x) 
the complete elliptic integral of the first kind. Compar- 
ison with the exact form of the real-space interaction 
potential (Fig.QJ indicates that the scaling form J3J is 
accurate also in lower LLs n > 1. Note that this ap- 
proximate expression is shown primarily for reasons of 
illustration here. In the calculations presented below, we 
use the exact form of the interaction potential, obtained 
by a Fourier transformation of Eq. (0) . Apart from small 
oscillations around the approximate form (black line in 
Fig-EJ, the potential exhibits a plateau of width 2Rc su- 
perimposed on the bare e 2 /er Coulomb potential, which 
is retrieved at large distances. The plateau is due to the 
ring-like form of the electronic wave functions in the n-th 
LL with a radius on the order of Rc- This model per- 
mits a common description of interacting electrons for 
all LLs. The Hamiltonian approach proposed by Murthy 
and Shanka»i2i is therefore not restricted to the lowest 
LL, but has universal validity, provided inter-LL mixing 
is negligible. 

For the following sections, we will choose a system of 
units with Ti = Ib = 1 to simplify the notation. 



III. HARTREE-FOCK TREATMENT OF THE 
ELECTRON-SOLID PHASES 

Whereas the quantum nature of the electron-liquid 
phases, which display the FQHE, is not captured in 
a mean-field treatment of the model, the Hartree-Fock 
approximation gives reliable estimates for the energy 
of states with a modulated electron density such as 
CDWsi&2iiiI In this section, we follow the lines of 
Refsi&iS. The Hartree-Fock Hamiltonian turns out to 
be linear in the projected density operators p(q) 



H 



HF 



1 

q 



(g)(p(-q))p(q) 



(5) 



with u™ b (q) = u„(q) - Ep w «(p) e 



(P*1y-Py1*)/( nB A). 

A detailed derivation of Eq. (J5J is given in the Appendix. 
Because the effective interaction potential in the n-th LL 
v n {q) is isotropic, the exchange potential is simply pro- 
portional to the real-space direct potential. For the inter- 
mediate LLs n = 1,2, and 3, the Fourier transformation 
may be calculated exactly without the help of the approx- 
imate scaling formula J3J|, which was used in a previous 
study^i A more detailed discussion of this approximation 
and its validity may be found in Sec. VIE. One obtains 
for the exchange potential in n = 1, 2, and 3 



Un=i(q) 
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where J\{x) is the first-order Bessel function. Using 



pc 



with the volume of the primitive unit cell of the triangu- 
lar lattice A pc — \/3A B /2, one obtains for the cohesive 
energy of the M-electron bubble phase 



Eg h {n;M,v) 



M 



(6) 



where G/ 7^ are the reciprocal lattice vectors of the 2D 
triangular lattice. 



B. Unidirectional CDWs - Stripe Phase 



with the modified Bessel functions Ij(x). It can be 
seen from Eq. Q that the cohesive energy E co h = 
(Hhf) I 'N e i is a functional of the order parameter A(q) = 
{p{q))/n B A 



j?CDW 
^coh 



(n; v) 



2v ^ " 
q 



(«)|A(q)| : 



It is the special form of the effective interaction potential 
in higher LLs (n > 1) that is responsible for the CDW 
formation: because of the form factor F n (q) , the effective 
interaction potential (J2J) vanishes at non-zero values of 
the wave vector qo(n). It is therefore energetically favor- 
able for the order parameter to have a maximum at these 
wave vectors, which leads to a density modulation in real 
space with a characteristic periodicity A ~ 2ir/q (n). 



A. Triangular CDWs - Bubble Phase 

We first investigate a CDW state of triangular sym- 
metry, which is also called bubble phase. This phase 
has been proposed as a candidate for the ground state 
at v < 1/2 in higher LLs. 8 ! 9 ' 10 Clusters of M elec- 
trons form a super Wigner crystal to minimize the resid- 
ual Coulomb interaction, and because these M-electron 
bubbles may be treated as classical objects, one ex- 
pects the CDW to have a triangular symmetry^ The 
phase is described by a local guiding-center filling fac- 
tor D(t) — Ylj ©(re — |r — R-jj), where r B = \/2M is 
the radius of a M-electron bubble and Rj are the lattice 
vectors of the 2D triangular lattice. Simple geometri- 
cal considerations lead to the relation r B /A B = \fH>v /2n 
with the lattice constant A^. Fourier transformation of 
the local filling factor yields the order parameter of the 
bubble phase 



A&(q) 



Aq 



-Ji(qV2M)^e lq - R3 , 



Around half-filling in higher LLs, a uni-directional 
CDW (stripe phase) competes with the bubble phase be- 
cause the latter would break the particle-hole symmetry, 
which is exact at v = 1/2. Other patterns such as a 
checkerboard pattern may also account for this symme- 
try, and there is a theoretical prediction that stripes be- 
come unstable towards the formation of an anisotropic 
Wigner crystal at T = 0. 23 Here, the discussion of the 
CDW phases is limited to the triangular and the uni- 
directional cases, which occur to be the most relevant 
phases observed in experiments. 

The stripe phase is described by a local guiding-center 



filling factor P(r) = V ■ ©(a/2 



), where a is the 



width of one stripe and Xj = jAs with integral j. The 
stripe periodicity As and the stripe width are obviously 
related to the partial filling factor v = a/ As- As for the 
bubble phase, the order parameter is obtained by Fourier 
transformation of v(r), 



AS(q)= 2^ o Sinfei^^ eiw . ASi 

^x Qx j 



where L x is the extension of the system in the ^-direction. 
This yields the cohesive energy for the stripe phase 



Eco h (n-A s ,9) 



n B 
2ir 2 D 



ly£0 



HP' 



27T 

q =AS l 



sin 2 (7rz^) 
P ' 



The optimal stripe periodicity is obtained from a mini- 
mization of the cohesive energy dE^ oh (n; As, v)/dAs = 
0. One then finds A s = 2.8AR C for n = 1, A<? = 2.76R C 
for n = 2, and As — 2.74Rq for n = 3. In the large-rt 
limit the optimal stripe periodicity, which is essentially 
independent of P, converges to the value As = 2.7 Rc in 
agreement with previous studies by Fogler et &1& 



IV. QUANTUM-LIQUID STATES AT 

P = l/(2s + 1) 

The incompressible quantum-liquid states at v = 
l/(2s + 1) in the lowest LL are described to very good 
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accuracy by Laughlin's wave functions^ A Laughlin- 
type state in an arbitrary LL may be obtained from the 
ansat»24 

vn! 

i 

where |f2 s ) is the Laughlin state in the lowest LL, and a\ 
is the usual ladder operator for the i-th particle connect- 
ing the different LLs. The energy of the Laughlin state 
in the lowest LL is usually defined with respect to the 
uncorrelated liquid and may be written as 

q 

where s(q) = (/>(— q)p(q))/N e i is the projected static 
structure factor. In order to obtain the energy of the 
corresponding state in an arbitrary LL n, it is suffi- 
cient to replace i>o(q) by the effective interaction po- 
tential v n (q) from Eq. @. The structure factor is 
the Fourier transform of the pair distribution function 
g s {r) = J z 3 ...z N (zi = 0,z 2 = r, z 3 , z N \tt s ), which is 
given by 21 i 25 

m=0 ' \ / 

The projected structure factor thus reads 

oo 

s(q) = (l-u)+4v^2 c^ m+1 L 2m+ i(q 2 )e-'? 2 / 2 ! 

m=0 

and the energy of the Laughlin state in an arbitrary LL 
n is 

U =-JxY, v ^ +E coh(n;s). (7) 
q 

The second term in Eq. (JjJ is the cohesive energy of the 
Laughlin state 

_ oo 

E h coh (n- s) = V - 4 m+1 V 2 n m+1 (8) 

m=0 

where the interaction potential J2J has been ex- 
panded in terms of Haldane's pseudopotentials, 
V? m+1 = (27rM)E q ^(g)Wi(9 2 )exp(-q 2 /2). 26 
The coefficients c| m+1 may be obtained from a fit of 
the pair-distribution function of the Laughlin states to 
Monte-Carlo calculationsASLSI An alternative, analyt- 
ical method to obtain these coefficients is to use the 
physical properties of the state described by Laughlin's 
wave functions, which result in a certain number of sum 
rules determining the coefficients 18 

CO 

/ j L 2m+1 — 2 ' 
m=0 



^(2m + 2)c| m+1 = --, 

oo 2 

^(2m + 3)(2m + 2)c| m+1 = i- ! 

m=0 

c s 2m+1 = -I for to < s. 

The first three sum rules are due to charge neutrality, 
perfect screening and compressibility, respectively. The 
last condition is given by the electron repulsion at short 
distances, where g(r) ~ r 2(2s+i) ^ ue ^ Q ^ ne j as trow fac- 
tors in the Laughlin wave function^ This method yields 
the coefficients 
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and C2 m+1 = for larger m. In spite of its simplic- 
ity, this method provides energy values for the quantum- 
liquid states which deviate less than 1% from Monte- 
Carlo calculations^ This accuracy is sufficient for the 
present comparison of the competing phases whose rela- 
tive energy difference is on the order of 10%. 

V. QU ASIPARTICLE /-HOLE EXCITATIONS AT 

P^l/(2s + l) 

Away from the filling factors vl — l/(2s + 1), the fi- 
nite energy of the excited quasiparticles (for v > vl) and 
quasiholes (for v < &l) must be taken into account. Al- 
though the quasiparticles/-holes carry an electric charge, 
they are treated in the following as non-interacting parti- 
cles. This approximation is valid for low quasiparticles/- 
hole densities, i.e. in the vicinity of vl = l/(2s + 1). One 
obtains for the cohesive energy of the quantum-liquid 
phase 

E q J(n; s, v) = E L coh {n- s) + [v{2s + 1) - l]A n (s), (9) 

where A n (s) is the energy of quasiparticles of charge 
l/(2s + 1) (quasiholes of charge — l/(2s + 1)) in units 
of the electron charge. 

The quasiparticle/-hole energies may be calculated an- 
alytically within the Hamiltonian theory proposed by 
Murthy and Shankar>iS The Hamiltonian is investi- 
gated in a CF basis, where the CF consists of an electron 
and a vortex-like excitation of vorticity 2s and charge 
— c 2 = —2ps/(2ps + 1), in units of the electron charge. 
If one neglects the internal structure of the CF, its elec- 
tric charge e* = (1 — c 2 )e leads to a reduced coupling 
to the external magnetic field. These CFs populate CF- 
LLs as do electrons in an external magnetic field. The 
integer p determines how many CF-LLs are completely 
filled, and the FQHE at fillings v — p/ (2ps + 1) may be 
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interpreted as an IQHE of CFs. For the present inves- 
tigations of quasiparticles/-holes around these fillings, it 
is sufficient to use the "preferred" combination for the 
density operators^ 

P P (q) = P(q) ~ c 2 x(q), 

where x(q) is the density operator of the vortex- like ex- 
citations ( "pseudovortex" ) . Because of its charge — c 2 , 
the pseudovortex density satisfies similar commutation 
relations as the original electron density operators ©, 



[%(q),x(k)] = -2isin 



(q x k); 

2c 2 



x(q + k), 



whereas both densities commute, [ / o(q),x(k)] = 0. At 
v = 1/2, which corresponds to the limit p — > oo and thus 
to c 2 = —1, one obtains the algebra proposed by Pasquier 
and HaldaneiSS 

The particular choice of the preferred combination re- 
spects the commutation relations @ in the small-wave- 
vector limit, whereas higher order corrections at larger 
wave vectors are strongly suppressed by the Gaussian 
in the effective interaction potential @. In terms of CF 
creation and annihilation operators ct m and Cj, m , respec- 
tively the density operator becomes^ 



p*(q) = n^ R \ni')(j\p p m')cl m c r 



with the matrix elements (m > m') 




iq y )l\ 



x U' 



V2 

2;*2 



q*l 



and (j > f) 



//! / -i{g x - iq y )l*B c 



ri-f I 1 'fl c 



„2(l-j+j' ) e -q 2 /2c 2 T j-j' ( q i*B 



L J 7 
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2c 2 



where l* B — 1 / V 1 — c 2 is the CF magnetic length and 
L™(x) are associated Laguerre polynomials. The quan- 
tum number j describes the CF-LL, and m is the CF 
guiding-center number. The ground state of the theory 
at v = p/(2ps + 1) consists of p completely filled LLs 
and is therefore characterized by the expectation value 

( c nj c j\m') = fij,j'fim,m'®(p— 1 ~~ j)- I n the Hartree-Fock 
approximation, the quasiparticle energies are thus given 
by the expression 

A qp(s,p) = (Cp, m fln4,m> - fin), 

and the quasihole energies are 

A qh( s >P) = { c l-i, m H n c P -i,m) - (H n ), 



where one averages over the ground state with the help 
of Wick contractions. This yields 



A" (s,p) 



«n(g)(p|^(-q)^(q)|p> 



r p_1 

/ v n ( q )j2\(p\p p m'W 

Jq j'=o 



(10) 



and 



A^(s,p) = ~Jv n ( q ){p-l\pP(-d)^(q)\p-l) 

+ I v n (q)J2\{P-nP P (<l)\3')\ 2 , (11) 
Jq i'=o 

where f = (l/-4)5Z q = J d 2 q/{2ir) 2 . The expressions 
(|10|l and ||TJ are generalizations to an arbitrary LL of 
Murthy and Shankar's results for n = 0~ In the lowest 
excited LLs n = 1, 2 and 3, one obtains the energies of the 
quasiparticle excitations for the Laughlin series (p = 1) 





s = 1 


s = 2 


s = 3 


s = 4 


n = 1 


0.2267 


0.1868 


0.1550 


0.1316 


n = 2 


0.1903 


0.1728 


0.1543 


0.1376 


n = 3 


0.1691 


0.1560 


0.1453 


0.1342 



and the energies of the quasihole excitations 



A^(s,P=l) 


s = 1 


s = 2 


s = 3 


s = 4 


n = 1 


-0.07172 


-0.07032 


-0.05887 


-0.04959 


n = 2 


-0.07876 


-0.07853 


-0.06728 


-0.05765 


rt = 3 


-0.07720 


-0.07944 


-0.07004 


-0.06124 



in units of e 2 /els- 



VI. THE DIFFERENT PHASES AND THE 
REENTRANCE OF THE IQHE 

Due to quasiparticle/-hole localization, the quantum- 
liquid phases display the FQHE at fillings ul = l/(2s + 
1), i.e. a plateau in the Hall resistance at values R xy = 
h/e 2 (N + v£) accompanied by a vanishing longitudinal 
resistance. On the other hand, the M-electron bubble 
phases are insulating because they are pinned by the un- 
derlying impurities. This suppresses the collective sliding 
mode, which, in the absence of a pinning potential, would 
result in a non-zero contribution to the electrical trans- 
port. The conduction is therefore entirely realized by the 
electrons in the completely filled lower LLs, and this leads 
to an integer-quantized Hall resistance R xy — h/e N. 
Whereas both the Wigner crystal and the bubble phases 
exhibit the IQHE, the Hall resistance for the stripe phase 
is not quantized and retrieves its classical value. 

We compare the cohesive energies of these different 
phases for the LLs n — 1,2, and 3. This allows one 
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FIG. 2: (a) Cohesive energies of the M-electron bubble, 
stripe, and quantum-liquid phases for n — 1 in units of e 2 /eZs- 
(b) RIQHE in upper spin branch of n — 1 measured by Eisen- 
stein et aJ^ The dot-dashed tangent on the curves M = 1 and 
M — 2 shows the energy of a mixed phase of 1- and 2-electron 
bubbles. 



to determine the ground state at different partial filling 
factors, and the results are compared to experimental 
findings. The discussion is limited to the range 9 < 1/2 
because the regime 9 > 1/2 is related to the former by 
particle-hole symmetry. The observation of this symme- 



try m experiment; 



,2.3.4 



supports the validity of the ap- 



proximation to treat the electrons in completely filled 
LLs as inert and the neglection of their spin degrees of 
freedom. 



A. The Phases in n = 1 

Fig-Hk shows the results for the cohesive energies in 
n = 1 in the absence of impurities. At low 9 < 0.23, the 
quantum-liquid phases are of a lower energy than the 
electron-solid phases. In the experiments by Eisenstein 
et al.£ the FQHE is observed around 9l = 1/5, whereas 
they find an insulating behavior below 9 < 0.15 corre- 
sponding to a magnetic field above B ~ 3.9T (Fig.JSJj). 
A 1-electron bubble is energetically favorable between 



0.23 < 9 < 0.30. Its insulating behavior is experimen- 
tally unveiled by the RIQHE around 3.8T. The quantum 
liquid has again a lower energy than the electron-solid 
phases between 0.30 < 9 < 0.36, and the FQHE is ob- 
served around 9l = 1/3. This alternation of electron- 
solid phases and FQHE states at Vl = 1/3 and 1/5 is 
in agreement with recent numerical studies in the frame- 
work of the density-matrix renormalization group^ The 
electron solid, however, appears in the numerical inves- 
tigations only in form of a strong oscillation of the pair 
correlation function and has been identified as a stripe 
phase. Our energy studies show that it is the bub- 
ble phase, which is energetically favored in this range 
of filling factors, and that the stripe phase would be 
the phase of lowest energy only at higher 9. In addi- 
tion, a mixed phase of 1- and 2-electron bubbles has a 
lower energy than the stripe phase up to a filling factor 
9 ~ 0.4 (dot-dashed tangent in FigJ^Ji). This tangent 
represents the convex envelope of the energy curves of 
the 1- and 2-electron bubble phase. Around 9 — 1/3 the 
quantum-liquid is still energetically more favorable than 
such a mixed phase. The experimental observation of the 
RIQHE around 3.62T supports the stability of a bubble 
phase at these fillings. Contrary to the presented theo- 
retical investigations, there is experimental evidence for 
a stripe phase around 9 = 1/2 only in the presence of 
an in-plane magnetic fields In the absence of such an 
in-plane field, which breaks the rotational symmetry of 
the system, a quantum-liquid phase is observed in form 
of the even-denominator FQHE states at v = 5/2 and 
7/2 in the n = l-LLAi^Si The origin of these enigmatic 
states is not completely understood until now. Theo- 
retical proposals range from spin-singlet formation^ to 
a Pfaffian wavefunction^ This wavefunction, which de- 
scribes a BCS-like state, may be a consequence of a pair- 
ing instability of the CF Fermi surface at 9 — 1/2^ A 
discussion of these states is, however, beyond the scope 
of the present article. 



B. The Phases in n = 2 

In the n — 2-LL, the energy studies (Fig.^Ji) show 
that the quantum-liquid phases are again favored at low 
densities 9 < 0.15 and between 0.19 < 9 < 0.21. For 
0.15 < 9 < 0.19, it is the 1-electron bubble phase which 
has the lowest energy. The 2-electron bubble phase is 
realized in a range 0.21 < 9 < 0.35. The stripe becomes 
the phase of lowest energy above 9 ~ 0.35, which leads 
to the strong anisotropy in the longitudinal magneto- 
resistance observed in experiments. 2 3 Fig.03 shows ex- 
perimental results obtained by Cooper et al& for the 
magneto- resistances in the lower spin branch of n = 2. 
Apart from the anisotropy around 9 = 1/2 (v = 9/2), 
a RIQHE is observed around B = 2.36T and 2.6T cor- 
responding to the filling factors 9 = 3/4 and 1/4. The 
presented energy calculations suggest that this insulating 
behavior is due to the formation of a 2-electron bubble 
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FIG. 3: (a) Cohesive energy for n = 2 in units of e 2 /e/_g. (b) 
RIQHE in lower spin branch of n = 2 observed by Cooper 
et al(- insets are a zoom on Hall (c) and longitudinal (d) 
resistance around B — 2.65T*- The dot-dashed tangent on 
the curves M = 1 and M = 2 shows the energy of a mixed 
phase of 1- and 2-electron bubbles. 



phase. The satellite maxima in the longitudinal resis- 
tance at B ~ 2.32T and 2.65T, accompanied by small 
minima in the Hall resistance (Figs. 01, d), indicate an 
incipient quantum melting towards the quantum-liquid 
phases around v — 1/5 or 1/7. Whether the 1/5 state 
is observable is not clear from the present energy studies 
because a mixed phase of 1- and 2-electron bubbles has 
approximately the same energy as the quantum liquid as 
is shown by the broken line in Fig.f5]}. 



C. The Phases in n = 3 

In the n = 3-LL, the quantum liquid is neither energet- 
ically favorable around P = 1/3 nor at v = 1/5 (Fig.@J. 
Also around v = 1/7 a mixture of 1- and 2-electron bub- 
bles may have a lower energy than the quantum-liquid 
phase, whereas a quantum melting may be observable at 
v = 1/9 in extremely pure samples. The stripe phase 
becomes stable above v = 0.36 in agreement with exper- 
imental observations^ 
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FIG. 4: Cohesive energy for n = 3 in units of e 2 /ds- The 
dot-dashed tangent shows the energy of a mixed phase of 1- 
and 2-electron bubbles, and the dashed tangent the energy of 
a mised phase of 2- and 3-electron bubbles. 



D. General Aspects of the Phase Diagram 

With increasing n the quantum-liquid phases shift to 
lower values of the partial filling factor, whereas the 
CDW phases become energetically favorable over a larger 
range of the electronic density in the last LL. This ef- 
fect may be understood from the scaling form of the ef- 
fective interaction potential in real space, as shown in 
Fig-H if the average distance d between electrons in the 
last LL is smaller than the width 2Rc of the plateau, 
it costs only a small amount of energy to decrease the 
distance between two electrons below d. At the same 
time a rather large energy on the order of the height of 
the plateau may be gained if the electrons cluster, thus 
reducing the number of other electrons with which they 
interact strongly^ The condition d < 2Rc leads, with the 
help of d ~ 1/ \f3 and Rq = \/2n + 1, to a critical value 
v c {n) oc l/(2n + 1) above which CDWs are expected to 
have a lower energy than the quantum- liquid phases. 16 
This scaling argument is supported by our energy inves- 
tigations, as well as by recent numerical studies based 
on the time-dependent Hartree-Fock approximation. In- 
deed, Cote et al~ found that in the n-th LL bubble 
phases with up to M = n + 1 electrons per bubble are 
energetically possible. Inspection of Figs.|5Ji, Et, and fj] 
confirm that as n increases, bubble phases with higher M 
are realized up to M = n+1. However, because large- M- 
bubble phases appear around v = 1/2, the bubble phase 
with n + 1 electrons is unstable towards the formation of 
a stripe pattern. 



E. Comparison with the Approximate Exchange 
Potential 

In a recent Rapid Communication^^ we used an ap- 
proximate form for the exchange potential in the calcu- 
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potential (dashed lines) and the results calculated with the 
exact exchange potential (solid lines), (a) results for n = 1; 
(b) results for n — 3. 



lations of the electron-solid phases. The expression 



4e 2 



-Re 



K 



1 - y/1 -4(2n+l)/V 



-i 2 



(12) 

is derived in the same manner as the scaling form of the 
real-space interaction potential Q and becomes exact 
in the large-n limit. For completeness, a comparison of 
the results in n = 1 and n — 3 obtained with the help 
of the exact and the approximate exchange potential is 
given in Fig.|SJ As expected, the approximation yields 
more accurate results in n — 3 than in n = 1, but it 
captures the correct physical properties also in the lower 
LLs. The rather compact form of the approximate ex- 
change potential 112(1 may therefore be an advantage in 
further theoretical analyses in intermediate LLs. 



VII. COMPETITION OF THE PHASES IN THE 
PRESENCE OF IMPURITIES 

Although the samples used for experimental investiga- 
tions of quantum Hall systems are nowadays extremely 



pure, the underlying residual impurities play an impor- 
tant role for the physical properties of these systems. On 
the one hand, they localize electrons/holes around v = N 
(IQHE) and quasiparticles/-holes around v — p/(2ps±l) 
(FQHE) and thus lead to the observed plateaus in the 
Hall resistances at these fillings. On the other hand, they 
pin the electron solids and therefore suppress their col- 
lective sliding mode. 

Whereas a weak impurity potential does not change on 
the average the cohesive energy of the quantum-liquid 
phases due to their incomprcssibility, the CDW phases 
can profit from the potential landscape by deformation of 
their charge structure. The overlap of the wave functions 
of electrons in different bubbles is negligible, and the bub- 
bles may therefore be treated as classical objects with 
an electric charge Me. The elasticity constant thus ap- 
proaches its classical value— fj, « 0.25M 2 e 2 n^ 2 /e, where 
n m = v /2-kM is the bubble density. 

We consider a short-range Gaussian impurity poten- 
tial with correlation length £ and strength V$. In the 
weak-pinning limit jak the energy gain for the M-electron 
bubble phase of elasticity [i may be obtained from a min- 
imization of the energy density 35 36 



m£_ 2 

L'n 



V, 



T 2 



(13) 



with respect to L . This yields the Larkin length 



L 



1 



47r eVo 



and the cohesive energy of the bubble phase © is thus 
lowered by the quantity 



6Eg h (M,V) 



e(L ) ^ (2tt) 3 / 2 Vi/e 



n e i 



fM&l* e 2 /e 



(14) 



The weak-pinning limit requires that Lq be larger than 
the lattice constant A^ , and the crossover to the strong- 
pinning regime may thus be characterized by the con- 
dition Lq < Ab- For the 1-electron bubble phase the 
strong-pinning case is equivalent to mere single-particle 
localization, and the energy gain per particle is simply 
given by SE coh = -V . 

The results for the cohesive energy of the bubble phases 
in n = 1 are shown in Fig.[S] for two different values of 
the impurity strength Vb/£ in comparison to the pure 
case. One observes that bubble crystals with larger num- 
bers M of electrons per bubble are less affected by the 
impurity potential due to their higher stiffness. Further- 
more, the quantum-liquid phases cease to be the energet- 
ically favored states in the dilute limit with v — > 0, and 
in agreement with experimental results^ we find that 
the Wigner crystal (1-electron bubble phase) is formed. 
The exact determination of this transition point, how- 
ever, requires a detailed knowledge of the nature and the 
strength of the underlying impurity potential. In addi- 
tion, correlation effects within the Wigner-crystal phase 
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FIG. 6: Cohesive energies of the 1- and 2-electron bubble 
and quantum-liquid phases for n = 1 in the presence of a 
short-range Gaussian impurity potential (in units of e 2 /els)- 
Full line: without impurities; dashed line: impurity strength 
Vo/£ = 0.005e 2 /el%; dashed-dotted line: impurity strength 
Vo/l = 0.010e 2 /e/|. 



which modify the energy curves should be taken into ac- 
count, as pointed out in previous studies for the lowest 



VIII. POSSIBLE EXPERIMENTAL 
INVESTIGATIONS 

In this section, we propose some new experimental in- 
vestigations, which could verify the presented theoretical 
results. The energy calculations indicate that the tran- 
sitions between the quantum-liquid and electron-solid 
phases are first-order quantum phase transitions. One 
consequence of this result is that the phenomenon of 
"super-cooling" is expected with respect to the control 
parameter, which is the filling factor tuned by the mag- 
netic field. Starting a transport measurement e.g. from 
the quantum-liquid phase at v = 3 + 1/5 and then lower- 
ing the magnetic field, the phase transition to the bubble 
phase would take place at a lower magnetic field than for 
the inverse case in which one raises the field starting from 
the bubble phase. The observation of such a hysteretical 
behavior in the Hall resistance may support the presence 
of first-order quantum phase transitions. 

Another possible investigation of the bubble phase is 
via transport studies under microwave irradiation of the 
2D electron system^ Due to the pinning of the bubble 
crystal by residual impurities, the longitudinal conduc- 
tivity 3?[(T xa; (w)] depends on the microwave frequency and 
exhibits a resonance at a finite value u p . The functional 
dependence of uo p was recently calculated for the Wigner 
crystal^ These results may be generalized to the bub- 
ble phase if one replaces the elasticity constant fi and the 
Larkin length Lq of the Wigner crystal by their values for 
the bubble phase derived in the preceding section. De- 
pending on the strength of the magnetic field, one finds 



for a classical bubble crystal 



„ „ M-V4p-3/4 5 for ^ fidds 

and 

^ „ _5° „ M -^ 2 D-^ 2 , for strong fields, 

where ITo ~ h/Lq is the characteristic self-energy of the 
phonons of the bubble crystal. The number of electrons 
per bubble may therefore in principle be determined from 
the resonance in the frequency-dependent conductivity 
at uj p . Because of the mixed phase with coexisting M- 
and (M + l)-bubbles, which is energetically favorable at 
certain filling factors, one expects to observe two peaks 
in $i[<r xx (uj)] whose respective weight varies with chang- 
ing magnetic field. This effect may be easier to measure 
around quarter-filling of the n = 3 than in n = 2 because 
in this range of v there is no competition with quantum- 
liquid or stripe phases. The resolution of such a double- 
peak structure, however, depends on the strength of the 
magnetic fields 

Recently, Cote et al. argued that the first- 
order quantum phase transitions between bubble crys- 
tals of different M would lead to experimentally ob- 
servable discontinuities in the magnetic susceptibility 
X = —{l/vnBA){d 2 E/dB 2 ) associated to the orbital 
magnetization^* It is, however, uncertain how these dis- 
continuities will evolve if the sharp phase transition is 
covered by an energetically favored mixed phase of bub- 
bles with different M. 



IX. CONCLUSIONS 

We have presented detailed energy calculations of 
the competing quantum- liquid and electron-solid phases, 
which allow for the determination of the ground state 
in the intermediate LLs n — 1,2, and 3 as a function 
of the partial filling v of the last LL. Whereas the ener- 
gies of the electron solids (triangular and uni-directional 
CDWs, i.e. M-electron bubbles and stripes), are cal- 
culated in the Hartree-Fock approximation, which gives 
reliable energy estimates for states with a modulated 
electron density^S**** the quantum nature of the liquid 
phases is not captured in a mean-field treatment of the 
model. 

At filling factors &l — V(2 S + 1), the quantum liq- 
uid may be described by Laughlin's trial wave functions. 
With the help of a set of sum rules imposed on these wave 
functionsj2ii2& we derived the structure factor of a quan- 
tum liquid, which has the same physical properties as 
the Laughlin liquid. 18 The determination of the structure 
factor allows one to compute the energy of the quantum- 
liquid phase, which agrees to great accuracy with numer- 
ical resultsi^i Away from precisely T>l = l/(2s + 1), the 
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energy of the quantum liquid increases due to the ex- 
citation of quasiparticles/-holes of finite energy A ra (s), 
which may be calculated analytically in the Hamiltonian 
theory of the FQHEiiS This non-monotonous behavior 
of the quantum-liquid energy gives rise to an alternation 
between quantum-liquid phases, which are energetically 
favored around vl = l/(2s + l), and electron-solid phases 
having a lower energy in between. 

Whereas the quantum liquids display the FQHE, the 
electron-solid phases, which form a triangular lattice 
of M-electron bubbles, are insulating because they are 
pinned by underlying impurities in the sample. The al- 
ternation of the different phases results in the observed 
RIQHE between the FQHE states at v = 1/5, 1/3, and 
the even-denominator state at half-filling in n = lX In 
n = 2, this competition between the different phases is 
manifest in the RIQHE found around v = 1/4 and 3/4 
and in the observation of an incipient quantum melt- 
ing around v ~ 1/7.. .1/5^ It is not clear whether the 
quantum melting may lead to a FQHE at v = 1/5 or 
1/7 in n = 2. Our energy calculations indicate that 
the 1/7 state is more stable than the 1/5 state for a 
pure sample. Impurities, however, favor electron-solid 
phases preferentially at lower v. The interpretation of 
the RIQHE in n = 2 as being due to the formation of an 
insulating bubble phase is further supported by recent 
microwave experiments, which reveal a peak in the lon- 
gitudinal conductivity at finite frequency^ This peak is 
likely to be caused by the excitation of a collective depin- 
ning mode>i£i2£iSi Microwave experiments in the RIQHE 
regime would thus clarify the origin of these insulating 
phases and could determine the number of electrons M 
per bubble because of the M-dependence of the resonance 
in the frequency-dependent conductivity. However, to 
the knowledge of the authors, such experimental investi- 
gations remain to be performed. 
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APPENDIX A: HARTREE-FOCK HAMILTONIAN 

In terms of fermionic operators, the Hamiltonian Q 
reads 



q ya,y' a 



X c 



- q x / 2 So + 9x / 2 C j/^ + <j x / 2 So ~ / 2 ' 



where the indices n at the operators are skipped for bet- 
ter legibility. This Hamiltonian becomes in the Hartree- 



Fock approximation 8 

Hhf = 5$>„(q) J2 e iqv{vo ~ v ' o) 
q yo,y' 

X [( C i ~q*/2 C yo+<lx/2} C y' +qj2 C y' -qJ 

~Ho-W2 c yQ-9*/2> c ^ +gx /2 c yo+gx/2 
= Hjj — Hp. 



The Hartree term is directly given in terms of the density 
operator 

Hh = i^u ff (q)(p(-q))p(q) 



with the Hartree interaction potential ujy(q) = v n (q). 
The average (p(q)) is proportional to the order parame- 
ter of the mean-field theory. In order to write the Fock 
term in the same manner, we perform the variable trans- 
formations 



Px = yo- y' , R = 



vo + y'o 



and in a second step 

J 2 ' N B 

y+,p y 



-E 



ipy(y+-y--q*) _ 



where N B — n B A is the total number of states per LL. 
Thus 



Hf = J5>n(q)$>^ 



(vo-y'o) 



;/(>•!/„ 



x ^ c \ -q I /2 c y[ ) -^/^) c \' a +q I /2 c yo+q^/ , 2 

nE^ E eWX ( C L 9x /2+ Px /2 C ^ 9 ,/2- Px /2) 



X C i?.+ 9x /2-p x /2 C «+9x/2+p x /2 



2N B 



Py.Px 



x ( c l_ +p x /2 c y- -p* /^ c l+ - Px /2 c y+ +Px /2 ■ 
We finally obtain 

Hf = ~^2u F (p){p(-p))p(p), 
p 

and one finds that the Fock potential is the Fourier trans- 
formed Hartree potential accompanied by an interchange 
of the x- and y-axes 



That both the direct and the exchange potentials are 
related to each other by a Fourier transformation, has 
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already been pointed out in previous works for the case of 
an isotropic interaction*^ However, if the interaction is 
anisotropic, the rotation of the frame of reference by ir/2 
has to be taken into account in the calculation of the Fock 
potential. The Hartree-Fock Hamiltonian can therefore 
be written entirely in terms of the order parameter and 
the projected density operator 

^ = 5$>ffF(q)<p(-qMq), (Al) 
q 



with the effective Hartree-Fock potential 

u hf (q) = u H (q) - u F {q) 

= «n(q) - -^5>„(p)e^-^).(A2) 
p 
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